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ABSTRACT 

Target  altitude  estimation  is  an  important  problem  in  over- 
the-horizon  radar  (OTHR)  systems  for  target  tracking  and 
classification.  Narrow  signal  bandwidth  and  rich  propa¬ 
gation  environment  in  OTHR  render  this  problem  difficult 
and  challenging.  In  this  paper,  we  develop  a  novel  technique 
that  allows  accurate  altitude  estimation  for  a  maneuvering 
target  during  its  ascending  or  descending  flight.  The  pro¬ 
posed  technique  considers  a  multiple-input  multiple-output 
(MIMO)  OTHR  platform  and  utilizes  the  difference  in  the 
Doppler  signatures  generated  in  a  micro-multipath  propa¬ 
gation  environment.  The  separability  of  the  multi-compo¬ 
nent  Doppler  signatures  in  the  time-frequency  domain  is 
used  to  individually  estimate  the  elevation  angles  of  these 
paths,  which  are  otherwise  infeasible  to  obtain,  due  to  their 
small  angular  separation.  The  estimated  elevation  angle  in¬ 
formation  is  then  used  to  achieve  a  reliable  estimate  of  the 
target  altitude. 

1.  INTRODUCTION 

Over-the-horizon  radar  (OTHR)  systems  are  critical  for  na¬ 
tional  defense  and  other  missions  because  of  their  wide-area 
surveillance  capabilities  to  cover  long  ranges  that  are  well 
beyond  the  limit  of  conventional  line-of-sight  (LOS)  radars 
[1]— [2] .  Skywave  OTHR  systems  exploit  the  reflective  and 
refractive  nature  of  high-frequency  (HF)  radiowave  propa¬ 
gation  through  the  ionosphere  to  cover  a  surveillance  area 
up  to  thousands  of  kilometers. 

Because  of  the  narrow  signal  bandwidth,  the  range  reso¬ 
lution  of  an  OTHR  system  is  usually  limited.  In  particular, 
target  altitude  information  is  difficult  to  obtain  and  remains 
a  challenging  problem  in  OTHR.  This  information,  even  in 
the  order  of  a  thousand  meters,  can  be  very  helpful  in  iden¬ 
tifying  the  type  of  aircraft  and  determining  the  flight  levels 
of  commercial  aircrafts.  Although  there  have  been  various 
attempts  to  estimate  the  absolute  or  relative  altitude  in¬ 
formation  [3,  4,  5],  absolute  target  altitude  estimation  in 
OTHR  has  not  received  proper  attention. 

A  practical  and  tangible  method  for  providing  the  al¬ 
titude  information  is  based  on  micro-multipath  model  that 
makes  use  of  multipath  returns  due  to  the  ocean  or  ground 
reflections  local  to  the  target  [3].  The  micro- multipath 
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propagation  is  illustrated  in  Fig.  1  with  a  simplified  flat- 
earth  ionosphere  model  [4].  Changes  in  the  target  altitude 
will  alter  the  distance  of  the  direct  path  and  the  multi- 
path,  generating  Doppler  variations,  with  opposite  signs, 
around  the  nominal  Doppler  signature.  These  propagation 
paths  of  the  emitted/received  signals  associate  themselves 
with  different  non-linear  time-frequency  trajectories,  each 
correspond  to  a  Doppler  signature  of  the  target  along  a 
propagation  path. 

It  is  shown  in  [4]  that  high-resolution  time-frequency 
analysis  of  the  Doppler  signatures  can  provide  resolution  of 
the  Doppler  signatures,  enabling  fast  and  robust  tracking  of 
the  altitude  changes.  However,  the  non-linearity  and  close 
separation  of  the  Doppler  structures  still  render  analyses  of 
such  signals  difficult  and  challenging  [4,  5].  Recent  works  [6, 
7]  have  developed  new  approaches  for  accurate  parametric 
estimations  of  time-varying  multi-component  signals  with 
closely  separated  Doppler  signatures  encountered  in  OTHR 
systems. 

In  this  paper,  we  utilize  successful  separation  of  the 
multi-component  Doppler  signatures  due  to  micro- multipath 
propagation  in  a  MIMO  OTHR  platform  to  estimate  the  ab¬ 
solute  altitude  of  a  maneuvering  target  during  its  ascending 
or  descending  flight.  It  is  noted  that  the  separability  of  the 
multiple  signal  components  in  the  time-frequency  domain 
has  several  advantages,  including  estimation  of  the  direc¬ 
tions  of  closely  spaced  targets  that  otherwise  are  difficult 
to  resolve  based  only  on  spatial  processing  [8]. 

The  following  notations  are  used  in  this  paper.  A  lower 
(upper)  case  bold  letter  denotes  a  vector  (matrix).  E[-]  rep¬ 
resents  the  statistical  mean  operation.  (•)*,  (-)T  and  (•)**"  re¬ 
spectively  denote  complex  conjugation,  transpose,  and  con¬ 
jugate  transpose  (Hermitian)  operations.  ln  expresses  the 
n  x  n  identity  matrix.  Diag(x)  denotes  a  diagonal  matrix 
using  the  elements  of  x  as  its  diagonal  elements.  In  addi¬ 
tion,  CNxM  denotes  the  complete  set  of  N  x  M  complex 
entries. 

2.  SIGNAL  MODEL 

2.1.  OTHR  Micro-Multipath  Model 

Consider  an  OTHR  system  using  the  flat-earth  model  as 
illustrated  in  Fig.  1.  In  practice,  the  transmit  and  receive 
array  are  separated  by  a  distance,  which  is  much  smaller 
than  the  range,  to  mitigate  power  leakage.  We  assume  that 
the  difference  of  the  range  due  to  the  array  positions  is  neg¬ 
ligible.  Thus,  the  slow-time  representation  of  the  received 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

JUN  2012 

2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2012  to  00-00-2012 

4.  TITLE  AND  SUBTITLE 

Altitude  Estimation  Of  Maneuvering  Targets  In  Mimo  Over-The  Horizon 
Radar 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory,AFRL/RYMD, Day  ton, OH, 45433 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS (ES) 

10.  SPONSOR/MONITOR’S 

ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

Presented  at  the  IEEE  Sensor  Array  and  Multichannel  Signal  Processing  Workshop,  Hoboken,  NJ,  June 
17-20,  2012 


14.  ABSTRACT 

Target  altitude  estimation  is  an  important  problem  in  overthe-horizon  radar  (OTHR)  systems  for  target 
tracking  and  classification.  Narrow  signal  bandwidth  and  rich  propagation  environment  in  OTHR  render 
this  problem  difficult  and  challenging.  In  this  paper,  we  develop  a  novel  technique  that  allows  accurate 
altitude  estimation  for  a  maneuvering  target  during  its  ascending  or  descending  flight.  The  proposed 
technique  considers  a  multiple-input  multiple-output  (MIMO)  OTHR  platform  and  utilizes  the  difference 
in  the  Doppler  signatures  generated  in  a  micro-multipath  propagation  environment.  The  separability  of  the 
multi-component  Doppler  signatures  in  the  time-frequency  domain  is  used  to  individually  estimate  the 
elevation  angles  of  these  paths,  which  are  otherwise  infeasible  to  obtain,  due  to  their  small  angular 
separation.  The  estimated  elevation  angle  information  is  then  used  to  achieve  a  reliable  estimate  of  the 
target  altitude. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

18.  NUMBER 

19a.  NAME  OF 

OF  ABSTRACT 

OF  PAGES 

RESPONSIBLE 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report 

(SAR) 

5 

PERSON 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


Consider  the  L  —  3  multipath  components  as  described 
earlier,  which  correspond  to  two  direction-of-departure 
(DOD)  angles  Ot,kt(i)  and  two  direction-of-arrival  (DOA) 
h  angles  Orjkr(i),  i  —  1?  *  *  • ,  L.  For  simplicity  and  without  loss 

of  generality,  we  only  consider  the  elevation  angle  0  which  is 
closely  related  to  the  target  altitude.  Then,  the  signal  data 
received  at  the  receive  array  corresponding  to  the  range  cell 
H  is  expressed  in  the  slow-time  domain  as  the  following  Nr  x  T 

complex  matrix 

X(t)  =  Arr(t)A?S  +  N(t),  (5) 


Figure  1:  Micro- multipath  propagation  in  OTHR. 


signal  is  expressed  as 

y(t)  crs7(t)  +  z(t)  m  (Tspe3^  +  z(t),  (1) 

where  t  is  the  slow-time  index,  a2  is  the  power  of  the  trans¬ 
mitted  signal  waveform,  y(t)  =  pe3^^  is  the  complex  re¬ 
flection  coefficient  of  the  target,  and  z{t)  is  noise. 

Denote  h  and  h  as  the  propagation  distances  of  the 
two  paths  observed  at  the  arrays,  respectively.  Then,  the 
combination  of  the  direct  path  and  the  multipath  in  both 
forward  and  return  links  yields  the  following  three  combi¬ 
nations  of  the  two-way  propagation:  path  I  (A  :  A),  path  II 
(l 2  M2),  path  III  (l  1  :  I2  and  I2  :  h).  Therefore,  the  received 
signal  consists  of  the  following  three  components: 

3  3 

y(t)  =  ^2  +  z{t)  =  aapie,4,i{t)  +  z(t),  (2) 

i= 1  i=  1 


where  the  phase  term  is  related  to  the  propagation  paths 
by  -27T fc(lt,kt(i)  +  lT,kr{i))/c,  where  fc  is  the  car¬ 

rier  frequency,  c  is  the  propagation  velocity,  and  h,kt(i)  and 
lr,kr(i)  respectively  denote  the  slant  range  of  the  forward 
and  return  links  associated  with  path  i. 

The  slant  ranges  l\  and  I2  can  be  expressed  in  terms 
of  the  ground  range  distance  R ,  the  ionosphere  altitude  if, 
and  the  target  altitude  h,  as 

h  =  (R2  +  (2H  -  h)2) 1/2  ,  h  =  (i?2  +  {2H  +  hf) 1/2  . 

(3) 

Equation  (3)  is  used  for  numerical  evaluations  of  the  slant 
ranges  in  this  paper.  Nevertheless,  for  illustrative  purpose, 
we  also  use  the  following  approximations,  which  hold  true 
when  R  H  ^$>  h,  to  reveal  a  clear  relationship  between 
the  slant  ranges  and  the  other  parameters: 


R  + 


2 if 2  -  2 iff 
R  ’ 


7  „  2if2  +  2Hh 

I2  ~R+ - ^ - •  (4) 


2.2.  MIMO  Radar  Signal  Model 

We  use  a  bistatic  MIMO  radar  model  for  the  direction 
finding  of  the  OTHR  system,  which  operates  in  a  pseudo¬ 
monostatic  mode.  The  MIMO  radar  consists  of  Nt  closely 
spaced  transmit  antennas  and  Nr  closely  spaced  receive 
antennas.  Denote  S  as  the  narrowband  waveform  matrix 
which  is  assumed  to  be  orthonormal  in  the  fast-time  do¬ 
main,  i.e.,  SS7/  —  I Nt- 


where  A*  =  [at(0t,fct(i)),  •  •  • ,  at(0t, **(£,))]  and  Ar  is# 
[ar(0r,fcr(i)),-  •  •  ,a r(0r,fer(L))],  with  at  (i) )  G  CNrXl  and 
eLr(Or,kr(i) )  G  CNtXl  respectively  denoting  the  transmit  and 
receive  steering  vectors  corresponding  to  the  zth  two-way 
path.  In  addition,  T(t)  —  Diag[7i(t),  •  •  •  ,7 z,(t)].  Moreover, 
we  assume  that  the  target  has  a  high  nominal  Doppler  fre¬ 
quency  that  clearly  distinguish  itself  from  the  clutter  region. 
The  signal  located  in  the  clutter  region  may  be  excluded 
for  direction  finding.  Thus,  N(f)  £  CNrXT  is  additive  noise 
matrix,  whose  elements  are  assumed  to  be  independent  and 
identically  distributed  (i.i.d.)  complex  Gaussian  random 
variables  with  zero  mean  and  variance  cr2.  To  justify  ex¬ 
pression  (5),  it  is  assumed  that  the  steering  vectors  remain 
unchanged  during  the  entire  slow-time  processing  period, 
which  is  often  the  case  for  far-held  targets. 

By  post-multiplying  (5)  by  SH,  we  obtain  Y(t)  £  xNr 
as 

Y(t)  =  Arr(i)A?  +  Z(t),  (6) 

where  Z (t)  =  N(t)SH. 


2.3.  Target  Elevation  Angle  and  Doppler  Signatures 

The  elevation  DOD/DOA  of  the  two  paths,  path  I  and  path 
II,  can  be  expressed  as 


0 1  =  tan 


2 H-h\ 
R  )’ 


62  =  tan 


(7) 


The  angular  separation  between  paths  I  and  II  is  usually 
very  small  and  simultaneous  estimation  of  these  angles  us¬ 
ing  high-resolution  direction  finding  techniques  would  thus 
be  infeasible. 

The  maneuvering  of  a  target,  in  general,  consists  of  hor¬ 
izontal  and  elevation  movements.  For  our  purpose,  we  con¬ 
sider  the  Doppler  frequency  characteristics  of  the  target  ma¬ 
neuvering  in  the  two  different  dimensions.  For  a  maneuver¬ 
ing  target,  we  consider  that  both  R  and  h  are  time- varying, 
whereas  the  height  of  ionosphere,  H ,  is  considered  constant 
over  the  processing  time. 

The  Doppler  frequencies  of  the  three  different  paths  are 
then  obtained  as 


,  ^  2/c  dh(t)  _  2/c  ^  4 fcH  M 

foAt)  =  ~  dT  ~  bWcv^ 

f°At)  =  ~~t - (i), 

fo,s(t)  =  fnAt)  =  k<m+!m  „  2JlK(t)vR(t), 

c  at  c 


(8) 


where  K(t)  =  1  —  2 H2 /R2 (t),  VR(t)  —dR(t)/dt  is  the  target 
velocity  in  the  range  direction  toward  the  radar,  and  vc(t)  = 
dh{t)/dt  is  the  ascending  velocity  of  the  target. 

From  the  above  discussion,  it  is  evident  that,  while 
the  dominant  Doppler  component  ]d  —  2 f CK (t)v R{t) / c  is 
shared  by  all  the  three  paths  and  reveals  the  target  velocity 
in  the  range  direction,  the  small  Doppler  difference  between 
the  paths,  fn  —  ±4 fcHvc(t)/[R(t)c],  is  a  function  of  h(t). 
The  Doppler  frequency  difference  allows  separation  of  the 
multi-component  Doppler  signatures  so  that  only  a  single 
component  is  selected  for  DOD  and  DOA  estimations. 

3.  TARGET  ALTITUDE  ESTIMATION 

Several  techniques  have  been  developed  to  discriminate  tar¬ 
gets/sources  based  on  time- frequency  signatures  (see  e.g., 
[6,  9]).  In  this  paper,  we  assume  that  perfect  Doppler  esti¬ 
mation  is  achieved  for  all  the  multipath  components. 

Rewrite  (6)  as 

L 

Y(t)  =  J2  V'Piit)  AieJ'*(t)  +  Z(t),  (9) 

i=  1 

where  A i  —  a-r(Orikr(i))&-i* Stationarizing  Y (t)  us¬ 
ing  the  estimated  phase  ^(t),  k  —  1,2,  which  is  assumed 
to  be  perfectly  estimated,  yields 

Y(t)e-^*(*>  =  aspk(t)  Afc 

+  +Z(t)e~jMt).  (10) 

i=l,i^k 

Therefore,  the  desired  signal  component  corresponding  to 
the  kth  path  becomes  a  stationary  DC  signal,  whereas  the 
other  multipath  components  will  have  non-zero  frequencies, 
provided  that  fu,i  ^  fD,k  for  i  ^  k.  As  such,  averaging 
the  above  result  over  the  observation  slow-time  amounts 
to  low-pass  filtering  the  signal,  which,  ideally,  will  capture 
the  desired  component  and  filter  out  the  multipath  compo¬ 
nents  represented  by  the  second  term  in  the  right-hand  side 
of  the  above  expression.  Further,  the  averaging  operation, 
given  below,  reduces  the  noise  power  proportionally  to  the 
coherent  integral  time.  That  is, 

1  To_1 

Y[fc]  =  Y(t)e_***(t)  =  aspk(t) Ak  +  Z[k] ,  (11) 

0  t= 0 

where  To  is  the  total  number  of  slow-time  samples  used  for 
the  operation,  and  =  (1/To)  o  1  Z (t)e~j^k^  is  the 
noise  component  after  filtering.  All  the  spatial  signature 
of  the  kth  path  is  carried  in  A k  in  the  above  expression. 
With  the  other  multipath  components  removed  or  highly 
mitigated,  and  the  noise  reduced,  accurate  DOD  and  DOA 
estimation  of  the  kth  path  becomes  possible.  This  proce¬ 
dure  is  repeated  for  k  —  1  and  2.  For  better  isolation  and 
filtering,  a  proper  window  (e.g.,  Hamming  window)  is  used 
to  reduce  the  sidelobe  leakage  of  undesired  signal  compo¬ 
nents. 

The  joint  DOD/DOA  estimation  in  a  bistatic  MIMO 
radar  system,  in  general,  may  involve  a  two-dimensional 
search.  A  number  of  techniques  have  been  developed  to 
simplify  the  two-dimensional  search  into  decoupled  one¬ 
dimensional  direction  finding  problem  and  perform  pairing 


of  the  DOD  and  DOA  information  when  multiple  targets 
are  present  (e.g.,  [10,  11,  12]).  In  the  underlying  problem, 
the  DOD/DOA  of  only  a  single  path  is  estimated  each  time. 
Therefore,  no  pairing  procedure  is  required.  In  this  paper, 
we  use  the  MUSIC-based  technique,  where  the  DOD  and 
DOA  are  estimated  from  Y^. 

Once  the  elevation  angle  of  paths  I  and  II  are  estimated 
as  6i  and  02,  we  can  readily  obtain  the  nominal  altitude  of 
the  target  from  Fig.  1  using 

h  =  i  [Z2  sin(02)  -  h  sin(0i)]  .  (12) 

The  above  operation  can  be  performed  for  both  DOD  and 
DOA  estimates.  Notice  that  the  actual  target  altitude  is 
time- varying,  and  the  above  estimate  reflects  its  mean  value 
during  the  observation  period.  The  instantaneous  target  al¬ 
titude  can  be  obtained  by  incorporating  its  relative  altitude 
changes  addressed  for  example  in  [4,  5]. 


4.  SIMULATION  RESULTS 


As  an  example,  we  consider  a  maneuvering  aircraft  which 
makes  a  360°  circular  turn  of  radius  5  km  in  approximately 
T  =  179.5  seconds  to  descend  by  approximately  2,250  me¬ 
ters.  The  centers  of  the  transmit  and  receive  arrays  are 
set  as  the  coordinate  origins.  The  other  key  simulation  pa¬ 
rameters  are  listed  in  Table  I.  All  the  multipath  signals  are 
considered  to  fall  within  the  same  range  cell.  The  aircraft 
maintains  a  horizontal  velocity  of  175  m/s  (630  km/hr).  Its 
elevation  velocity  varies  sinusoidally,  and  the  corresponding 
target  altitude  is  expressed  as 


h{t)  =  h{  0) 


VC, 


,max  T0 


7T 


(13) 


Table  1:  Key  Parameters 


Parameter 

Notation 

Value 

initial  range 

R(  o) 

1,500  km 

ionosphere  height 

H 

160  km 

aircraft  initial  height 

MO) 

10,000  m 

horizonal  target  velocity 

VR,  max 

175  m/s 

maximum  descending  velocity 

Vc,  max 

19.68  m/s 

carrier  frequency 

fc 

16  MHz 

waveform  repetition  frequency 

fs 

40  Hz 

The  actual  Doppler  signatures  obtained  from  (3)  are  de¬ 
picted  in  Fig.  2(a),  and  those  obtained  from  the  data,  after 
the  path  I  component  is  made  stationary  using  its  exact 
Doppler  signature,  are  shown  in  Fig.  2(b).  The  maximum 
value  of  the  one-side  Doppler  difference  is  0.448  Hz.  The 
first  and  last  20  seconds  are  excluded  from  the  DOD/DOA 
estimation  because  of  the  close  Doppler  signatures  between 
the  different  multipath  components.  A  Hamming  window  of 
length  5601,  which  spans  the  entire  selected  time  period  of 
approximately  140  seconds,  is  used  to  mitigate  the  sidelobe 
leakage. 

For  simplicity  and  without  loss  of  generality,  we  con¬ 
sider  small-size  arrays  that  use  6  transmit  antennas  and  10 
receive  antennas.  Both  arrays  are  linear  and  the  antennas 
are  extended  in  the  range  direction.  To  achieve  a  large  ar¬ 
ray  aperture,  the  minimum  redundancy  array  (MRA)  con¬ 
figurations  [13]  are  used  at  both  the  transmit  and  receive 
arrays,  where  the  unit  separation  is  one  wavelength. 


In  Fig.  3,  the  MUSIC  pseudo  spatial  spectra  are  shown 
for  10  independent  estimates  respectively  obtained  for  DOD 
and  DOA,  where  the  input  SNR  is  —10  dB.  In  Fig.  4,  the 
root  mean  square  error  (RMSE)  of  the  mean  target  alti¬ 
tude  versus  the  input  SNR  is  obtained  by  averaging  over 
100  independent  trials.  It  is  evident  from  these  results 
that  the  elevations  angles  of  the  multipath  components  are 
clearly  resolved  for  target  altitude  estimation.  The  DOA 
estimates  are  generally  of  better  quality  when  compared 
with  the  DOD  estimates  because  more  sensors  are  used  at 
the  receiver.  The  RMSE  of  the  mean  altitude  estimate  is 
improved  by  an  order  of  magnitude  as  the  input  SNR  is 
increased  by  20  dB.  Note  that  a  1  km  RMSE  is  achieved  for 
the  small  array  configurations  even  when  the  input  SNR  is 
below  —20  dB. 


(a)  DOD  (b)  DOA 


Figure  3:  DOD  and  DOA  estimates. 


5.  CONCLUSION 

We  have  developed  a  novel  technique  for  the  estimation  of 
a  maneuvering  target  altitude  in  a  MIMO  OTHR  system. 
By  utilizing  the  difference  between  the  multi-component 
Doppler  signatures  due  to  micro-multipaths,  the  proposed 
technique  is  based  on  individual  estimation  of  the  DOD 
and  DOA  information.  Further,  the  knowledge  of  the  non¬ 
linear  Doppler  signatures,  which  can  be  estimated  by  us¬ 
ing  recently  developed  time-frequency  analysis  techniques, 
enables  them  to  become  stationary  for  achieving  SNR  en¬ 
hancements  through  extended  coherent  integration  times. 
Simulation  results  have  shown  that  the  proposed  technique 
is  very  promising  in  a  realistic  low  SNR  environment. 


Figure  4:  RMSE  of  the  target  altitude  estimates  for  differ¬ 
ent  input  SNR. 
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